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Abstract 

A "natural" model for the QCD invariant (running) coupling, free of the IR singularity, is 
0^ ' proposed. It is based upon the hypothesis of finite gluon mass m s \ existence and, technically, 

uses an accurate treating of threshold behavior of Feynman diagram contribution. The 
model correlates with the unitarity condition. 
, ^ ■ Quantitative estimate, performed in the one-loop approximation yields a reasonable 

lower bound for this mass m g i > 150 MeV and a smooth IR freezing at the level a s (Q 2 ) < 1. 

S: 

cnJ : 1 Introduction 

The issue of the infrared (IR) behavior of the strong interaction becomes more and more 
7— i ! actual from the physical point of view along with the further experimental data accumu- 

lation. In the perturbative quantum chromodynamics (pQCD) this behavior is burdened 
ff-j ' with "unphysical" singularities marked with the so called scale parameter A ~ 300 MeV. 

These singularities contradict some general principles of the local QTF. In the "small mo- 
. mentum transfer region" Q < 3A they violate the weak coupling regime and complicate 

theoretical interpretation of data. 
Qh! Quite recently attempts have been made to devise models for invariant (running) 

coupling a s (Q 2 ) and for observables free of singularities in the IR region. One idea was 
to exploit the scheme arbitrariness of the third beta-function coefficient (3$. Performing 
an "optimization" of the beta-function the authors of paper |l[ argued for special solution 
with an IR fixed point and a s (Q 2 ) < 1 bounded in the whole region Q 2 < lGeV 2 . 
In the paper || the allowance was made for nonperturbative contributions into a s {Q 2 ) 
emerged from the vacuum QCD background fields. Here, the adjusting parameter, the 
background (hybrid) gluon mass Mb, enters as an IR regulator InQ 2 — ^ ln(Q 2 + M|) of 
the gluonic logs. At Mb = 1.5 GeV the running coupling maximum value turned out to 
be a s (0) ~ 0.4 . One more trick || uses an imperative of the Kallen-Lehmann analyticity 
in the complex Q 2 plane. Effectively, it results in the smooth freezing of a s (Q 2 ) at the 
level of 0.5 - 0.7 for 0.1 GeV 2 < Q 2 < 1 GeV 2 and in the universal, i.e., loop-independent 
IR limiting value a™ ax = a s (Q) = 4tt /ft ~ 1.4. 

In the last, so called "invariant analytic", approach (for a recent review, see the 
finiteness and smoothness of the a s (Q 2 ) behavior in the IR region is achieved only due to 
the analyticity imperative, without introducing any adjustable parameters. Qualitatively, 
the effect of smooth freezing arises here due to the addition of a power (in the Q 2 variable) 
and nonanalytic (in the coupling constant a s ), i.e., nonperturbative terms which restore 
the Kallen-Lehmann analyticity for renormalization - invariant quantities. 
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In this note, we consider one more possibility of constructing the invariant coupling 
ots{Q 2 ) that is free of unphysical singularities and does not involve explicit nonperturbative 
contributions. Here, the Q 2 power contributions appear due to threshold effects and an 
essential physical ingredient is the assumption of the finite gluon mass existence. 

Our model expression for a s (Q 2 ) is obtained by the renormalization group (RG) 
summation of mass-dependent one-loop Feynman diagrams contribution — see, below, 
Eqs. (0) and fl3|). It depends upon the gluonic m gl and light quark m u d s masses, in the 
IR region Q 2 > obeys a nonsingular behavior with a finite limiting value a s (0) , and as 
Q 2 /m 2 — > oo smoothly transits into the usual asymptotic freedom formula. 

2 Massive loops 

Our starting point is quite simple and natural — we propose to take into the account the 
threshold mass dependence while considering the infrared region. As it is well known, the 
"leading UV logs" , which after the RG summation yield, in particular, the Landau pole, 
arise from the one-loop Feynman diagrams. However, in the IR region these diagrams 
behavior, generally, is far from being logarithmic. For instance, to the virtual dissociation 
of a vector particle (photon, gluon) into a massive fermion-antifermion pair (e + + e~ ; q + 
q) in the s-wave state, there corresponds a function I S (Q 2 /M 2 ) that can be represented 
in the form of a spectral integral 

OO 

1 

and in the space-like region z > is a positive, monotonically growing function with 
logarithmic asymptotic behavior I s (z) ~ Inz — C s + 0(z~ 1 ) ; C s = 5/3. For definiteness 
and in accordance with the QED tradition we have subtracted it at z = . This agreement 
turns out to be convenient for our purpose in the "low Q 2 region" . 

In QCD, besides I s (z) , essential is a function I p (z) describing a virtual dissociation 
of a vector gluon into a pair of massive vector or scalar particles (gluons or "ghosts" ) in 
the p-wave state. This function can also be represented in the form (p]) with an adequate 
weight function k p (a) and obeys the same simple properties at z > . In particular, in 
the UV limit I p (z) ~ Inz — C p + 0(1/ z). In our analysis only asymptotic constant C p will 
be of importance. It relates to the integral over the phase volume. Due to this C p > C s . 
For the qualitative estimate we use C p = 8/3. 

Here, we assume that a virtual gluon has a finite (reasonably small — see below) mass 
m g i . We postpone for the future any detailed discussion of a possible origin of this mass 
noting that the mechanism should be based upon a deeper understanding of the ground 
state structure of the quantum gauge SU(3) field. As a provisional ad hoc working model 
one could imply the picture of spontaneous symmetry breaking, analogous to the SU(2) 
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case (for a fresh discussion of the subject, see, e.g., Ref.|J). In addition, one should have 
in mind the upper phenomenological bound 

m g! < 600 MeV , 

corresponding to the absence of a direct experimental signal for this mass existence. 



3 Massive Renorm— group 

For an analysis of the invariant QCD coupling at small space-like values of the variable 
q2 < ^2 we uge " m assive", that is mass dependent, renormalization group as it has 
been explicitly formulated in the pioneer papers J7| in the mid-fifties. 

In particular, we shall exploit the fact that massive RG, quite in parallel to the widely 
used massless one, sums all iterations of a one-loop contribution in the invariant coupling 

a s (Q\ert = a s - a 2 s Ai(Q 2 , m 2 ) + a 3 s (A^Q 2 , m 2 )) 2 + . . . 
into the geometric progression M 

Here, we use the notation for the so-called couplant 

a s {Q 2 ) = a s {Q 2 )/A7r , a s = a s /Aix . 

Analogous (approximate) RG summed mass-dependent expressions are known M for 
the two-loop case — see below Eqs. @ and @ — as well. 

Let us make a comment on the renormalization (subtraction) scheme (RS). We use the 
subtraction at Q 2 = , that is the MOM-scheme instead of the massless MS one which 
is in common practice. Among QCD practitioners there exists a strong prejudice against 
any MOM-schemes due to their gauge dependence. For our analysis it is essential to use 
mass-dependent expressions for the diagram contributions. The particular RS is not of 
principal importance. One could transit from our scheme to a MS scheme by standard 
rules. In particular, connection between MOM, massive MS and the popular massless 



MS schemes has been discussed in detail in papers || [UJ and we have no possibility of 
repeating the discusson in this short note. 



4 The one-loop analysis 

For the one-loop contribution to an invariant couplant we use the expression^ 
MQ\m*) = n I. g ) - § £ g) . (3) 

1 Generally, the A\{Q 2 ,m 2 ) should combine several (propagator and vertex) contributions. However, 
for our semiquantitive analysis only UV behaviour with asymptotic constants C s ,p will be essential. 
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This expression by convention turns to zero at Q 2 = . Hence, in (§) a s = a s (0) . In the 
opposite limit at Q 2 ^> m 2 , we have 



In the three-quark region, this gives 



m 2 m 2 3 m u md 



O 2 

A!(Q 2 ,m 2 ) ~ 9 In -\ -2 In - - In 



26 ; Q 2 > m 2 , m; 



In the massless case to the expression l/a+v4i(Q 2 , m 2 ) there corresponds 9 ln(Q 2 /A 2 
Hence, we get the relation 



1 , m= A 
— + 22 In — + 18 In — 



a. 



m Rl 



4 , m 2 
26 + - In s — 

3 m u m d 



between the combination l/a— 221nm g i and the Q CD scale parameter A. For a quantita- 
tive estimate let us define an "effective one-loop QCD scale parameter" from the condition 
a s (M 2 ) = 0.37 that yields A : = 250 MeV. 

Now putting m u = 5 MeV, m d = 10 MeV, m s = 150 MeV , m g i = m s /a , we arrive at 
the relation 



— = 25 

a. 



22 In a 



a s (0) 



2tt 



12.5 + llln(m g i/ 



(4) 



The most important qualitative result that follows from it consists in the existence of 
a reasonably small lower bound for the gluon mass corresponding to the a s (0) ~ 1 
condition. 
We have 

100 MeV < m g i < 600 MeV . 
In the Table a few values of the IR limit a s (0) for some m g i of this interval are given. 

Table 



m gl /MeV 


100 


150 


200 


300 


450 


a.(0) 


0.87 


0.55 


0.41 


0.36 


0.31 



As it follows from the Table, there exists a rather wide interval of the gluonic mass 
m g i with "reasonably small" a s (0) < 1 values. 

Let us note that the lower bound as obtained from relation ([|) has only a qualitative 
nature. The point is that due to the relation m Uj d <C m g \ the derivative da(x)/dx at 
x = is positive and in the region < Q 2 < m 2 ^ the light quark contribution dominates. 
The real place of possible "blowing up" is close to Q 2 = 2m 2 j . As a result, the real lower 
bound for m g i corresponds to the condition a s < 1/ (28 — 4 In a) , and turns out to be close 
to the strange quark mass. 
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The upper bound existence for the a s (0) value resembles us the property of unitary 
models for lower partial waves of hadron scattering in the low-energy elastic region. These 



models were popular in the sixties - see, e.g., Refs. |jnj| . Besides analyticity, they satisfy 
the two-particle unitarity condition. One- loop diagrams summed in the solution @, (|3]) 
just relate to this last condition . This is the reason for associating our construction with 
unitarity. 



5 Two— loop corrections 

For a more accurate numerical description of the invariant coupling in the IR region one 
can use the two-loop massive perturbation expansion 

a s (Q 2 ) P ert,2 = a s - alA^Q 2 , m 2 ) + +a 3 s Aj - a 3 s A 2 (Q 2 , m 2 ) + . . . . (5) 

An approximate two-loop massive RG solution is of the form 

a s (Q 2 )r S ,2 = a s jl + a s Ai(Q\ m 2 ) + ^[^'^ ^ (l + ^(Q 2 , m 2 )) J . (6) 

At small a s values this expression corresponds to (||). At the same time, at Q 2 3> m 2 it 
can be represented in the usual form 

a7 l (Q 2 )r g , 2 - ft {in g + h In (in } ; b x = | . 

As it follows from (|6|), in the IR region 

= - + A^Q 2 , m 2 ) + a s A2(Q 2 , m 2 ) . (7) 
a s {Q ) a s 

Hence, the two-loop contribution A 2 (Q 2 ,m 2 ) is here suppressed by a small numeri- 
cal factor a s = a s (0)/47r < 1/10 and cannot seriously influence the one-loop estimate 
obtained above. 



6 Conclusion 

The estimate obtained shows that an accurate description of the threshold effects related 
to the light quarks and, especially, to gluons allows us to formulate one more resolution of 
the issue of unphysical IR singularities in the pQCD. The price of this resolving consists 
in introducing of the only parameter - the gluon finite mass. 

Thence, it is possible to keep the QCD invariant coupling a s (Q 2 ) in the weak coupling 
domain in the whole space-like region < Q 2 < oo for the gluon mass values 

200 MeV < m g i < 600 MeV . 
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This theoretical "window" for the gluonic mass could be enlarged if we change the light 
quark masses from their current values to some effective ones m e J[ ~ m n \TA\ . Here, one 
can use the gluon effective mass of the same order of magnitude. 

To conclude, let us remark that the present model, in reality, is not very far from 
our previous construction Refs. |3j with an explicit introducing of nonperturbative contri- 
bution. A gluon mass m g \, being considered in the light of the genuine gauge-invariant 
QCD Lagrangian, certainly represents an effective nonperturbative parameter. The same 
is true for the abovementioned rrfj^. 

In our opinion, this demonstrates once more (compare, e.g., with discussion of the 
QED case in Ref.Jl5|) that the ghost-pole trouble is not a physical problem. It is a 
technical drawback inherent to usual, the pQCD one, way of theoretical analysis. 
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